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We study the effects of jet quenching on the hydrodynamical evolution of the quark-gluon plasma 
(QGP) fluid created in a heavy-ion collision. In jet quenching, a hard QCD parton, before frag¬ 
menting into a jet of hadrons, deposits a fraction of its energy in the medium, leading to suppressed 
production of high-pT hadrons. Assuming that the deposited energy quickly thermalizes, we sim¬ 
ulate the subsequent hydrodynamic evolution of the QGP fluid. For partons moving at supersonic 
speed, Vp>Cs, and sufficiently large energy loss, a shock wave forms leading to conical flow Q. The 
PHENIX Gollaboration recently suggested that observed structures in the azimuthal angle distri¬ 
bution might be caused by conical flow. We show here that, for phenomenologically acceptable 
values of parton energy loss, conical flow effects are too weak to explain these structures. 
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Recent Au-|-Au collision experiments at the Relativis¬ 
tic Heavy Ion Collider (RHIC) saw a dramatic suppres¬ 
sion of hadrons with high transverse momenta (“high- 
Pt suppression”) and the quenching of jets in the 
direction opposite to a high-pT trigger particle 
when compared with p-l-p and d-|-Au collisions. This is 
taken as evidence for the creation of a very dense, color 
opaque medium of deconfined quarks and gluons 1^. In¬ 
dependent evidence for the creation of dense, thermal- 
ized quark-gluon matter, yielding comparable estimates 
for its initial energy density ((e) ^ lOGeV/fm^ at time 
'Ttherm ^ 0.6fm/c |3|), comes from the observation of 
strong elliptic flow in non-central Au-|-Au collisions 
consistent with ideal fluid dynamical behaviour of the 
bulk of the matter produced in these collisions. 

These two observations raise the question what hap¬ 
pens, in the small fraction of collision events where a hard 
scattering produces a pair of high-pr partons, to the en¬ 
ergy lost by the parton travelling through the medium. 
The STAR Collaboration has shown that, while in central 
Au-|-Au collisions there are no hard particles left in the 
direction opposite to a high-pr trigger particle, one sees 
enhanced production (compared to p+p) of soft (low-pr) 
particles, broadly distributed over the hemisphere dia¬ 
metrically opposite to the trigger particle . This shows 
that the energy of the fast parton originally emitted in 
the direction opposite to the trigger particle is not lost, 
but severely degraded by interactions with the medium. 
As the impact parameter of the collisions decreases, the 
average momentum of the particles emitted opposite to 
the trigger particle approaches the mean value associ¬ 
ated with all soft hadrons, i.e. the (pr) of the ther- 
malized medium 0. This suggests that the energy lost 
by the fast parton has been largely thermalized. Nev¬ 
ertheless, this energy is deposited locally along the fast 
parton’s trajectory, leading to local energy density inho¬ 
mogeneities which, if thermalized, should in turn evolve 
hydrodynamically. This would modify the usual hydro¬ 


dynamic expansion of the collision fireball as observed in 
the overwhelming number of soft collision events where 
no high-pT partons are created. 

Since the fast parton moves at supersonic speed, it was 
suggested in Ref. 0 that a Mach shock (“sonic boom”) 
should develop, resulting in conical flow and preferred 
particle emission at a specific angle away from the di¬ 
rection of the fast parton which lost its energy. This 
Mach angle is sensitive to the medium’s speed of sound 
Cs and thus offers the possibility to measure one of its key 
properties. A recent analysis by the PHENIX Collabora¬ 
tion 0 of azimuthal di-hadron correlations in 200 A GeV 
Au-I-Au collisions revealed structures in the angular dis¬ 
tribution which might be suggestive of conical flow. 

The idea of Mach shock waves travelling through com¬ 
pressed nuclear matter was first advocated 30 years ago 
111 ^ I but RHIC collisions for the first time exhibit |3 
the kind of ideal fluid behaviour which might make an 
extraction of the speed of sound conceivable. An alter¬ 
nate scenario, in which the color wake field generated by 
the fast colored parton travelling through a quark-gluon 
plasma accelerates soft colored plasma particles in the 
direction perpendicular to the wake front min , leads 
to an emission pattern which is sensitive to the prop¬ 
agation of plasma rather than sound waves na . In a 
strongly coupled plasma with overdamped plasma oscil¬ 
lations, which seems to be the preferred interpretation of 
RHIC data [HIH , the wake field scenario should reduce 
to the hydrodynamic Mach cone picture. We here study 
the dynamical consequences of the latter, going beyond 
the discussion of linearized h;^rodynamic equations in a 
static background offered in . 

We assume that just before hydrodynamics become ap¬ 
plicable, a pair of high-pT partons is produced near the 
surface of the fireball. One of them moves outward and 
escapes, forming the trigger jet, while the other enters 
into the fireball along, say, the —x direction. The fireball 
is expanding and cooling. The ingoing parton travels at 


2 


the speed of light and loses energy in the fireball which 
thermalizes and acts as a source of energy and momen¬ 
tum for the fireball medium. We model this medium 
as an ideal fluid with vanishing net baryon density. Its 
dynamics is controlled by the energy-momentum conser¬ 
vation equations 


= r, 


( 1 ) 


where the energy-momentum tensor has the ideal fluid 
form = {e+p)u^u'' with energy density e and 

pressure p being related by the equation of state (EOS) 
p = p{s), = 7 ( 1 , Vx, Vy, 0) is the fluid 4-velocity, and the 

source current is given by 


r{x) = Jix) {1,-1,0,0), 
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Massless partons have light-like 4-momentum, so the cur¬ 
rent describing the 4-momentum lost and deposited 
in the medium by the fast parton is taken to be light¬ 
like, too. rjet(t) is the trajectory of the jet moving with 
speed \dxjet/dt \ =c. ^{x) is the energy loss rate of the 
parton as it moves through the liquid. It depends on the 
fluid’s local rest frame particle density. Taking guidance 
from the phenomenological analysis of parton energy loss 
observed in Au-|-Au collisions at RHIC we take 


dE s{x) dE 
dx So dx 


( 4 ) 


where s{x) is the local entropy density without the jet. 
The measured suppression of high-p^ particle production 
in Au-|-Au collisions at RHIC was shown to be consistent 
with a parton energy loss of ^ |q = 14 GeV/fm at a ref¬ 
erence entropy density of sq = 140 fm“^ For compar¬ 
ison, we also perform simulations with ten times larger 
energy loss, = 140GeV/fm. 

For the hydrodynamic evolution we use AZHYDRO 
nm, the only publicly available relativistic hydrody¬ 
namic code for anisotropic transverse expansion. This 
algorithm is formulated in {T,x,y,g) coordinates, where 

T=^/t‘^—z'^ is the longitudinal proper time, In 

is space-time rapidity, and r_L = {x, y) defines the plane 
transverse to the beam direction z. AZHYDRO employs 
longitudinal boost invariance along z but this is violated 
by the source term ©• We therefore modify the latter 
by replacing the d-function in @ by 


^3(r-rjet(t)) 
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with cr = 0.35 fm. Intuitively, this replaces the “needle” 
(jet) pushing through the medium at one point by a 


“knife” cutting the medium along its entire length along 
the beam direction. The resulting “wedge flow” is ex¬ 
pected to leave a stronger signal in the azimuthal parti¬ 
cle distribution dN/dcf) than “conical flow” induced by a 
single parton, since in the latter case one performs an im¬ 
plicit (/)-average when summing over all directions of the 
cone normal vector. While a complete study of this would 
require a full (3-|-l)-dimensional hydrodynamic calcula¬ 
tion, the present boost-invariant simulation should give 
a robust upper limit for the expected angular signa¬ 
tures. We show that the angular structures predicted 
from wedge flow are too weak to explain the experimen¬ 
tally observed 0 -modulation Q- 

The modified hydrodynamic equations in (r, x, y, p) co¬ 
ordinates read Hfl 


+ dxiVxfn + dyiiyfn 

drf^^ + dx{Vxr^) + dyiVyf^^) 
drf^y + dxivxf^^) + dyivyf^y) 


-p + J, ( 6 ) 

-dxp - J, (7) 

-dyp, ( 8 ) 


where f^^'' = p = rp, and J = tJ. 

To simulate central Au-|-Au collisions at RHIC, we use 
the standard initialization described in ^ and provided 
in the downloaded AZHYDRO input file Ig , correspond¬ 
ing to a peak initial energy density of Eq = 30 GeV/fm^ 
at To = 0.6fm/c. We use the equation of state EOS-Q de¬ 
scribed in incorporating a first order phase tran¬ 

sition and hadronic chemical freeze-out at a critical tem¬ 
perature Tc = 164 MeV. The hadronic sector of EOS-Q is 
soft with a squared speed of sound « 0.15. 

In our study the quenching jet starts from xjet = 6.4 fm 
at tq =0.6fm, moving left towards the center with con¬ 
stant speed Ujet = c. For an upper limit on conical flow 
effects, the fast parton is assumed to have sufficient ini¬ 
tial energy to emerge on the other side of the fireball. 
To simulate cases where the fast parton has insufficient 
energy to fully traverse the medium we have also done 
simulations where the parton loses all its energy within 
an (arbitrarily chosen) distance of 6.4 fm. We further 
compared with a run where the parton moved at (con¬ 
stant) subsonic speed (ujet = 0.2 c< Cs). 

The resulting evolution of the energy density of the 
QGP fluid is shown in Figs.n]and|21 The left column of 
Fig.m shows results for a phenomenologically acceptable 
value ^ I p =14 GeV/fm [13 for the reference parton en¬ 
ergy loss whereas in all other columns we use a ten times 
larger energy loss. The width of the Gaussian source 
(see Eq. 0) is cr=0.7 fm. In the left column of Fig. Q] 
the effects of the energy deposition from the fast parton 
are hardly visible. Only for a much larger energy loss 
(right column) we recognize a clear conical flow pattern. 
The accumulating wave fronts from the expanding energy 
density waves build up a “sonic boom” shock front which 
creates a Mach cone. The right columns in Figs. Hand[3 
show that the cone normal vector forms an angle 0m with 
the direction of the quenching jet that is qualitatively 
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FIG. 1: Contours of constant local energy density in the x-y 
plane at three different times, t = 4.6, 8.6, and 12.6 fm/c. In 
each case the position of the fast parton, along with the in¬ 
tegrated energy loss AE = J J{x)dxdydT up to this point, 
is indicated at the top of the figure. Diagrams (a)-(c) in 
the left column were calculated with a reference energy loss 
di5/da:|o = 14 GeV/fm, those in the right column (panels (d)- 
(f)) with a 10 times larger value. 


consistent with expectations from the theoretical rela¬ 
tion COS0M = Cs/vjet- However, this angle is not sharply 
defined since the cone surface curves due to inhomogene¬ 
ity and radial expansion of the underlying medium. This 
differs from the static homogeneous case Q. 

When the parton travels at a subsonic speed 
Ujet = 0.2 c < Cs (left column in Fig. [21), it doesn’t get very 
far before the fireball freezes out due to longitudinal ex¬ 
pansion. In this case one only observes an accumulation 
of energy around the parton but no evidence of Mach 
cone formation. When the parton travels with njet = c 
but looses all its energy after 6.4 fm before fully travers¬ 
ing the fireball (right column in Fig. the fireball evo¬ 
lution beyond r=7 fm is not affected by the fast parton 
directly but only indirectly through the propagation of 
earlier deposited energy. Still, Figs.^ and|2f show that 
the late time evolution of the hreball is quite similar in 
both cases, demonstrating that energy deposition by the 


FIG. 2: As in FigQbut for a parton moving at subsonic speed 
■Cjet = 0.2 c (left column) and for a fast parton (njet = c) which 
loses all its energy within the first 6.4fm (right column). 

fast parton during the late fireball stages is small, due to 
dilution of the matter, and can almost be neglected. 

Tests with different values for the width cr of the Gaus¬ 
sian profile in Eq. for the deposited energy show that 
the cone angle gets better defined for smaller source size 
cr. Note that the quenching jet destroys the azimuthal 
symmetry of the initial energy density distribution but 
leaves the azimuthally symmetric energy contours to the 
left of the jet unaffected. 

We close with a discussion of observable consequences 
of conical flow. One expects 0 azimuthally anisotropic 
particle emission, peaking at angles (p = 7rzt 6 m relative 
to the trigger jet where 9m is the Mach angle. Using 
the standard Cooper-Frye prescription, we have com¬ 
puted the angular distribution of directly emitted pions 
at a freeze-out temperature Tfo = 100MeV 0. Figure (21 
shows the azimuthal distribution of 7r“ from a variety 
of different simulations. [The (/>-independent constant 
/dyd4>)-ao jet = 27 from central collisions without 
jets has been subtracted.] For dE/dx\Q = 14GeV/fm the 
azimuthal modulation is very small. In none of the cases 
studied we find peaks at the predicted Mach angle with 
an associated dip in the direction of the quenched jet at 
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FIG. 3: Azimuthal distribution dN/dydcp of negative pions 
per unit rapidity. In the upper panel we integrate over all pr 
while the lower panel shows only pions with pr > 1 GeV/c. 
Different symbols refer to different parameters as indicated. 
For better visibility the ^-independent rate in the absence of 
the quenching jet has been subtracted. Filled symbols show 
the realistic case dE/dx\o = 14 GeV/fm, enhanced by a factor 
10 for visibility. 

(f> = TT. One rather sees a peak at 4> = n, broadened by 
shoulders on both sides which turn into small peaks at 
(relative to the quenching jet) backward angles if thermal 
smearing is reduced by considering only high-pj’ pions. 
The peak at 0 = it is absent when the parton loses all its 
energy halfway through the fireball or is too slow to get 
to the other side before freeze-out, suggesting that it re¬ 
flects the directed momentum imparted on the medium 
by the fast parton. It is slightly more accentuated for 
smaller cr and higher pr- We also found that the width 
of the shoulders is almost independent of the speed of 
sound of the medium and can not be used to diagnose 
the stiffness of its equation of state. The shoulders exist 
even for subsonic parton propagation (ujet = 0.2 c, upright 
triangles), showing that other mechanisms (such as back- 
splash from the hard parton hitting the fluid and a gen¬ 
eral bias of the energy deposition towards the right side 
of the fireball due to the higher density of the medium at 
early times) have a strong influence on the angular dis¬ 
tribution of the emitted particles which interferes with 
the position of the Mach peaks. The shoulder resulting 
from this combination of effects is much broader than the 
angular structures seen in the data 0. The absence of a 
clear dip at ^ = ir in our simulations is all the more trou¬ 
bling since it should have been stronger for the “wedge 


flow” studied here than for real “conical” flow. 

Our calculation does not average over the initial pro¬ 
duction points of the trigger particle, i.e. it ignores that 
in most cases its quenching partner does not travel right 
through the middle of the fireball cylinder, but traverses 
it semi-tangentially. This should further decrease the 
prominence of the shoulders in dN/dcf). We conclude that 
conical flow may be able to explain the broadening of the 
away-side peak in the hadron angular correlation func¬ 
tion around (j) = TT pointed out by the STAR Collabora¬ 
tion Q , but is unlikely to be responsible for the relatively 
sharp structures near (f) = iril seen by PHENIX Q . This 
conclusion extends to other conical flow phe nomena, such 
as those generated by color wake fields [ll|,|l3 . 
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